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ABSTRACT 

A new mechanism of particle acceleration to ultra high energies, driven by the rotational 
slow down of a pulsar (Crab pulsar, for example), is explored. The rotation, through the time 
dependent centrifugal force, can very efficiently excite unstable Langmuir waves in the e-p plasma 
of the star magnetosphere via a parametric process. These waves, then, Landau damp on electrons 
accelerating them in the process. The net transfer of energy is optimal when the wave growth 
and the Landau damping times are comparable and are both very short compared to the star 
rotation time. We show, by detailed calculations, that these are precisely the conditions for the 
parameters of the Crab pulsar. This highly efficient route for energy transfer allows the electrons 
in the primary beam to be catapulted to multiple TeV (~ 100 TeV) and even PeV energy domain. 
It is expected that the proposed mechanism may, partially, unravel the puzzle of the origin of 
ultra high energy cosmic ray electrons. 



Subject headings: Particle acceleration 

1. Introduction 

If one could devise a mechanism that could 
convert even a tiny fraction of the enormous 
rotational energy of pulsars (spinning neutron 
stars) into particle kinetic energy, an outstand- 
ing problem of high energy astrophysics may 
get closer to a solution. One could, then, pro- 
vide, for instance, a possible explanation for 
what accelerates cosmic ray electrons to ul- 
tra high energies- a regime recently explored 
by th e H.E.S.S., Pamela and Fermi collabora- 
tions jAharonian et at.ll2008t IChang et at.l l2008t 
lAbdo et al.ll2009t lAckermann et al.ll2010[ ) 



Some of the most significant observations and 
developments in the realm of high energy cosmic 
ray electrons may be summarized as: 1) the 2008 
announcement of the H.E.S.S. collaboration show- 
ing evidence for a substantial steepening in the en- 
ergy s pectrum of cosm i c ray electrons above 600 
GeV ( Aharonian et at. 20081 ). 2) Suggestions for 
the possible sources for the very high energ y elec- 



trons (VHE electrons) IChang et at.l (|2008l ): it is 



conjectured, for example, that either pulsars or 
the annihilation of dark matter particles might 
be the source of ultra relativistic energies, 3) the 
launching of the Fermi spacecraft into a near-earth 
orbit on 11 June 2008 - the new instruments on 
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board detected cosmic-ray electrons up to ITcV 
and c onfirmed an excess of VHE leptons in cosmic - 



rays (|Abdo et all l2009t lAckermann et al.l l2010f) . 



and 4) Analysis of the observational data of the 
Fermi telescope on millisecond pulsars (the sur- 
face magnetic field close to B st ~ 10 8 5 G) by 
Kisaka fe Kawanaka ( 2012t ) positing that the data 
may have evidence of extremely high energy (~ 
50TcV ) cosmic ray electrons. 

This paper is an attempt to formulate a theo- 
retical framework for particle acceleration driven 
by the rotational slow down of a pulsar, in par- 
ticular, the Crab pulsar (B at - 7 x 10 12 G). The 
star rotational energy is channeled to the parti- 
cles in a two step process - the excitation of Lang- 
muir waves (via a two-stream instability) in the 
electron-positron plasma in the pulsar magneto- 
sphere; and the damping of Langmuir waves on 
the fast electrons accelerating them to even higher 
energies up to and beyond 100 TeV. 

All known pulsars are observed with decreasing 
spinning rates. The Crab pulsar, for example, ro- 
tates with a period P ps 0.0332s with a character- 
istic rate of change P = dP/dt ps 4.21 x 10~ 13 ss _1 . 
The slowing down releases rotational energy that 
must transform into some other form. The slow- 
down luminosity of the star is given by W f=s 
where I = 2MR 2 t /5 is the moment of inertia of 
the pulsar, R st ps 10 6 cm is its radius, fi = 2ir/P 
is its angular velocity, f2 = dfl/dt, and the pulsar 
mass M ~ 1.5 x M Q where M « 2 x 10 33 g mea- 
sures the solar mass. Out of this enormous output, 
W w 5.5 x 10 38 erg/s, only a part ~ 2 x 10 38 erg/s 
is sufficient to power the Crab Nebula; the rest - 
almost 60% of the energy budget - is released via, 
yet, unknown channels. 



In the standard pulsa r model (jDeutshl 11955 



Goldreich fc Julian Il969h . an electrostatic field 



uproots particles from the pulsar's surface. These 
particles are then accelerated along the magnetic 
field lines by the nonzero longitudinal electric field. 
The electrons, accelerated to relativistic energies, 
begin to emit copious curvature radiation. The 
emitted photons with energies beyond the pair 
production threshold 2mc 2 (m is the electron's 
mass and c is the speed of light), in turn, cre- 
ate electron-positron pairs thought the channel, 
7 + B —> e + + e~ + 7' (B is the dipolar magnetic 
field). Newly produced pairs are further acceler- 
ated and emit photons; The cascade lasts till the 



nascent pair plasma, especially in regions farthest 
from the star, is dense enough to screen out the 
initia l electrostatic field ( Sturrokl 1971 : Tademarul 



19731) 



This height h of the gap region (where the elec- 
tric field is nonzero) is, therefore, strictly lim- 
ited. The limited accelerating region, in turn, 
imposes an upper limit on the maximum attain- 
able electron energy. The potential difference 
in the gap has been estimated to be AV ps 



1.6 x W 2 B^ /7 P- 1 / 7 ^ (jRuderman fc Sutherland! 



19751 ). where B 12 = B/W 12 G 



6.7 is the di- 
mensionlcss magnetic induction. Thus the po- 
tential difference, predicted by the standard gap 
model, may boost the electron Lorentz factors 
to ~ 4 x 10 6 - not sufficient to explain the ob- 
served very high energy (VHE) radiation. Figure 
1 shows the schematic representation of the dis- 
tribution function of plasma particles in standard 
models of the pulsar magnetosphere. The primary 
beam particles are shown by the narrow shaped 
area with the Lorentz factor 7^. The wider region 
represents the distribution of secondary electron- 
positron pairs, characterized by a corresponding 
Lorentz factor, 7 P when the distribution function 
has a maximum value. 

Several mechanisms have been proposed to, 
somehow, enlarge the gap and boost up en- 
ergy accumulat io n in this zone. In particular, 
Usov fc Shabadl (|1985T ) proposed the so-called 



inte rmediate positronium form ation mechanism, 
and Arons fc Sharleman ( 19791 ) studied the field 
lines supposedly curved towards the rotation and 
rectifying und er certain conditions. Invok ing gen- 
eral relativity, iMuslimon fc Tsvean ( 1992 ) consid- 



ered the role of Kerr metric in the creation of the 
additional electric field. In all these attempts, the 
gap size does increase but not enough to account 
for the observations. 

Equally noteworthy are the attempts to under- 
stand high energy acceleration in terms of the 
Fermi acceleration mechan ism ( Blandfordl 1980t 
Blandford fc Ostrikenll980( ). The Fermi-type ac- 



celeration is found to be efficient when the seed 
population of pre-accelerate d electrons possesses 
quite high Lorentz factors (jRieger fc Mannheim 



2000h . The problem of reaching such high Lorentz 



factors, still, remains unresolved. 

It is in this backdrop - searching for a mech- 
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anism for particle acceleration- that we turn to 
rapidly rotating pulsars. The challenge will lie in 
uncovering a "new" mechanism for pumping out 
and harnessing energy from the central engine. 

We explore, here, a highly efficient mecha- 
nism for transferring energy from pulsar rotation 
to particle kinetic energy. The invoked mech- 
anism, first suggested (in a limited context) in 
(jMachabeli et all 120051) . consists of two princi- 
pal steps: 1) The time dependent centrifugal 
force, imparted to the relativistic plasma particles 
from star rotation, parametrically, drives Lang- 
muir waves, and 2) These waves damp on very en- 
ergetic electrons to accelerate them even to higher 
energies. In order to appreciate the workings of 
the mechanism, the following short review of the 
standard picture of the composition of the pulsar 
magnetosphere may be necessary. 

When acting on relativistic particles, the cen- 
trifugal force induces a periodic motion in time 



(jMachabeli et al.l l2005j ). and, as a result, ends 



up fuelling a parametric instability that may ef- 
ficiently pump energy from rotation into Lang- 
muir waves. The possibility of parametric reso- 
na nce in an electron plasma was fi rs t discussed 
bv IMachabeli k Tsikarishvilil (|l978l) : iMachabeli 
(|1978l ). It is worth noting that centrifugal ef- 
fects strongly influence particle acceleration and 
MHD processes i n the magnetospheres of pulsars 
(lOsmanov et al.l 120091: lOsmanov fc Riegerl liooj ) 
and AGN (jOsmanov et al.ll2007t )~ 

We will begin, here, by a radical and far- 
reaching gene ralization of the mecha nism, first 

" (|2005l ). The re- 



introduced in Machabeli et al 



fashioned model, applied to the Crab pulsar, 
departs from the very simplified approach of 



Machabeli et al.l (120051 ) that examined only a 
narrow subclass of centrifugally accelerated rel- 
ativistic electro ns- the ones with radial ve locity 
v (t) « c cos fit (jMachabeli fc Rogavalll994l ). We, 
instead, extend the initial particle velocity profile 
by letting Vo(t) ps ccos (fit + </>) with 0^0, the 
general solution allowed by the centrifugal force. 
It is somewhat unexpected and rather amazing 
that insisting on a nonzero initial phase spawns a 
rather strong instability. 

After this long introduction, we work out, in 
Section 2, the theory of centrifugally induced 
Langmuir waves. In section 3, we apply the model 
to the Crab pulsar and in section 4, we summarize 



our results. 

2. Centrifugally Driven Langmuir Waves 
in the Pulsar magnetosphere 

To develop a model for the centrifugally in- 
duced parametric Langmuir instability, we will 
assume that the magnetic field lines are almost 
rectilinear- a good approximation for distances less 
than the curvature radius of the field lines. Then 
in the rotating frame of straight field lines co- 
rotating with the pulsar, the metric is given by 



ds 2 



1 



fl 2 r 2 



(it 2 



dr l . 



(1) 



In the 1 + 1 formalism (jThorn et al.l Il986f ) for a 
single particle, the particle equation of motion, in 
the rotating frame of reference, reads 



dPi e, / 1 „ 

-T 1 = 7ig + - E + - Vi x B 

ar m \ c 



(2) 



where i = e,p is the species index, e$ is the 
charge of the corresponding particle, dr = £dt, 

1/2 

£ = (l — f2 2 r 2 /c 2 ) is the so-called time lapse 
function, g = — V£/£ is the gravitational poten- 

— 1/2 

tial, 7i = (l — V 2 /c 2 ) is the Lorentz factor, 
and pi and Vj are, respectively, the dimensionless 
momentum (p, — > pi/m) and velocity. The first 
term on the righthand side of the equation repre- 
sents the centrifugal force, and plays a crucial role 
in our model. 

To transport the equation of motion into the 
laboratory frame, we use the well known relation 
7 = £7' connecting Lorentz factors in the rotat- 
ing (7) and laboratory (7') frames. The identity 
d/dr = d/d(£dt) + (v • V) takes Eq. © into the 
laboratory frame 



dpi 

dt 



-c 2 7^V£ 



+ (vi • V)p 4 



E+ -v ; x B 

c 



(3) 



where we have omitted the primes for notational 
simplicity. Note that the first term, originating in 
the relativistic effects of co-rotation, has a singu- 
lar behavior on the light cylinder surface (a hypo- 
thetical zone, where the linear velocity of rotation 
exactly equals the speed of light). Throughout 
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the area of interest, the overall effect of the term 
c 2 7i£V£ will be significant. 

Equation <j3j> is clearly fluid-like. In order to 
motivate this change in appearance, we need to 
make a digression. The plasma we plan to de- 
scribe is really a kinetic fluid, the distribution 
function covers a wide range of 7 (energy) and 
as will later insist, a variable phase. A proper 
detailed description, then, is rather involved and 
will require numerical modeling. In this paper, 
however, we plan to extract our main results in a 
highly reduced description. We will assume that 
the plasma is multi-stream, each stream with a 
characteristic 7 and a characteristic phase. Then 
we will solve the linear "interacting" dynamics of 
two such streams and demonstrate how a very in- 
teresting manifestation of the two stream instabil- 
ity takes place; the phase difference between the 
streams will emerge as the driver. Each stream, of 
course, is an electron-positron plasma. Equation 
([3]) , then, is to be viewed as equation of motion of 
the specie i(electron or positron) in a given stream, 
and along with the continuity equation 



^ i + V-(n J v i ) = 
at 



and the Poisson equation 



V • E 



(4) 



(5) 



completes the dynamics. 

Th e centrifugal force (jMachabeli fc Rogava 
is the cause of equilibrium or background 
velocity for the stream particles. The elec- 
tromagnetic terms do not contribute (to the 
equilibrium) because of the frozen-in condition 
E n + vni x Bn = 0. In equi librium, Eq. d3]) reads 



(|Machabeli fc RogavallT99l 



a r 



1 



n 2 r 2 



1 - 



2 
72 



(6) 



For highly relativistic particles (7 ^ I), Eq. 
yields the simple periodic solution, 

r(t)«^sin(fit + 0), 

v (t) « V cos (Clt + 0) . 

The phase c/>, necessary to insure a general solu- 
tion of Eq. ([5]) , was taken to be zero in previous 



(7) 

(8) 



work IjMachabeli et al.ll2005l ). This was a very se- 
rious shortcoming because the zero phase solution 
applies only to a very small subset of particles. 
The nonzero phase will insure the inclusion of all 
particles in the dynamics. Consequently, the in- 
stability could feed on much of the population of 
the magnetospheric plasma resulting in very large 
growth rates . 

The orbits (|7l8p . along with the frozen-in con- 
dition, provide the equilibrium input for the lin- 
ear stability analysis. Since the calculated equilib- 
rium is time dependent, the standard time -Fourier 
anaylsis will not pertain- the solution of the sta- 
bility problem will be more involved and nuanced. 

One begins with the standard splitting of the 
physical variables = (n, v, p, E,) into equilib- 
rium and perturbed parts 



(9) 



The next step is somewhat different; the linear 
perturbations are expanded as 

^(t.r) oc^^exp^kr)], (10) 

where both the amplitude and the exponent are 
functions of time. To keep the notation simple, 
we will drop the superscript I labelling perturbed 
quantities. 

In the linearized version of Eqs. (J5][S|): 

^ + ikvopi = vott 2 rpi + —E, (II) 
at m 



drij 
dt 



+ ikviorii, +ikn i0 v 2 = 0, (12) 



ikE = 4tt^ 



eprii, 



(13) 



E is the electrostatic field, and i is the species 
index. Notice that we have retained only the radial 
component of the linearized momentum equation, 
Pi = e r ■ Pi and k = e r ■ k. 

The electron-positron dynamics can be com- 
bined in terms of two fluid variables:, the relative 
velocity v = v e — v p , the differential density?! = 
n e — n p , the average fluid velocity v° « v° = v°, 
and the mean density, n® = n® = n°. Thus for the 
two stream model (labelled by /3=1,2), the linear 
system (|1HI13[) takes the form, 



Spa ., „9 2e 

— — + ikvpopp = vpoll rpfj H E, 

at m 



(14) 
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—- + ikv /soup, +lkn/3oV / 3 



ikE 



(15) 
(16) 



the two streams being coupled through the Pois- 
son equation. We remind the reader that the two- 
stream model was invoked, very crudely, to cap- 
ture the essence of the vast spectrum in energy 
and initial orbit phase available to the plasma par- 
ticles. Each stream carries its own 7 and phase <p. 

Equations (|14j|16j) , due to the time dependence 
of the orbits (r, and Vq), represent a set of non- 
autonomous differential equations. The goal is to 
determine if an arbitrary perturbation will grow in 
time. We will follow two complementary method- 
ologies to establish the existence of an instability 
and to determine the growth characteristics; an 
exact numerical solution and an approximate an- 
alytical estimate via a quasi-modal approach. 

We convert (|14j|16p into two coupled equations 



d 2 N : 



dt 2 



2 + uj 2 2 N 2 = -uj 2 2 N ie - tx 



in terms of the perturbed "densities' 



(17) 
(18) 

(19) 



where we have replaced Vq by c. The frequen- 
cies wi,2 = J &ire 2 ni t2 /mjf 2 are > respectively the 
effective relativistic plasma frequencies (for longi- 
tudinal modes) corresponding to the two streams, 
and 



ck 



X 



fi 

with 2<j)± 



sin ( VLt + (f> 1 ) - sin ( M + (f> 2 ) ] 



± 



and 

, ck , 
= 2— sin( 



bcos(nt+(f> + ) 
(20) 

(21) 



In deriving Eqs. (|17H18|) . we have invoked (ck/Vl) » 
1 to drop the first term on the right hand side of 
Eq. (fT4|) . This parameter will turn out to be ap- 
proximately the ratio of the Langmuir frequency 
to the star rotation frequency and will lie in the 
range 10 6 - 10 10 . 

The evolution Eqs. (|17fll8() with definitions ([20} 
l2Tjl have the following features: 



1) The evolution of the two streams is coupled. 

2) The time dependence of the coefficients is 
fully contained in x- 

3) x depends on both the sum and difference of 
the phases associated with the two streams. 

4) The dependence on the phase sum 2<f> + = 
4>i + 4>2 is trivial. In fact it can be absorbed in a 
redefinition of time variable; <f> + will be dropped 
from now on. 

5) The phase difference 0_, however, is a funda- 
mental parameter of the two-strem system and de- 
termines the magnitude of x- For <j)_ = 0, x ~ 0, 
Eqs. (|17Hr5|) reduce to an autonomous system with 
constant coefficients and lead to the modal disper- 
sion relation uj 2 = uji 2 + uj 2 2 . The two streams col- 
lapse to one since (/>_, the "measure" of their dif- 
ference, has gone to zero. Strictly speaking, this is 
only an approximation, since in Eq. (|20[) we have 

assumed Vb/3 = c, whereas Vq/3 ~ c ^1 — • 
Without this approximation, we will get a nonzero 
growth rate even for 0_ =0). However, this ves- 
tigial growth rate is many orders of magnitude 
smaller than what we will obtain for finite 0_ . 

We would reassert that introducing the ini- 
tial phase is a fundamental departure from earlier 
work. It is the phase difference that differentiates 
the two streams and sets the stage for a strong two 
stream instability. 

We will first perform a time Fourier analysis 
to search for a growing quasi mode (for non au- 
tonomous systems there are no fourier modes). We 
must manipulate Eqs. tjTTl [TS)) to yield 



d 2 N 2 



U) 2 



,d 2 N 2 



dt 2 CJ! 2 dt 2 ' 

We, then, use the Bessel identity 

e ±i6 S ini, J„(6)e ±< «' ) 



(22) 



(23) 



and take the Fourier transform of Eqs. (|17H22j) to 
derive the relation 

_ (uj + (-1-1 + v)Slf 



(u 2 -co 2 1 )N 1 (co)=uj 2 2 J2QtQ 



fit; 



(24) 



where 



Mb) 



Equation (f2"4")l is clearly not a dispersion rela- 
tion, each Fourier component (labeled by uj) is 
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connected to all other components. In spite of the 
formidable double sum on the right hand side, we 
can still extract useful information that will estab- 
lish the possibility of the growth of perturbations. 
An indicative " dispersion relation "relation may 
be obtained by keeping terms with v = /i 

o, 2 -u,?-u, a a Jg(6) = a,' ]T 4(b) , f )2 . (25) 

We make further analytical progress by simply 
keeping a single resonant term defined by the har- 
monic of the star rotation frequency that is in res- 
onance with the real part of the mode frequency, 
u r = Hres^ where u = uj r + A ( A << u> r ). After 
some straightforward algebra, one can derive 

A 3 = ^ J ^ (i)2 , (26) 
implying an imaginary part- growth rate: 

The preceding analytical calculation, though sim- 
ple and revealing, should be taken only as an in- 
dicator that the Langmuir waves (charge separa- 
tion perturbations) can grow in time transferring 
energy from pulsar rotation to the electric field. 
As remarked earlier the instability is an appropri- 
ate version of the two-stream instability in which 
both streams move with speeds close to c, but 
whose constituent particles start with different ini- 
tial orbits-the difference manifesting itself through 
the phase factor <j>_. 

The expression (|27|) for the growth rate tells us 
that, apart from its dependence on uj\ and W2, the 
growth rate is controlled by the Bessel function 
J Prtl (i)) = J Mre< ((2cfc/fi)sin^)_) of very high or- 
der fj, res = w/fi 10 4 . Bessel functions of such 
high order arc vanishingly small unless the argu- 
ment (2ck/Q,) sincj)—) m /i res . For 0_ approaching 
zero, the growth rate becomes vanishingly small 
(jMachabeli et al.ll2005l) . 

We could readily make analytical estimates for 
the growth rate of the rotationally driven Lang- 
muir waves for a typical pulsar. However we must 
bear in mind that these estimates can be merely 
suggestive since the fundamental dynamics is non- 
autonomous in time and a given Fourier mode can- 
not capture the time behavior. It is particularly 



so because /i res is so large that a whole band of 
harmonics will contribute to the sum (f!H|) even 
if mode coupling were neglected. However, the 
evolution of each mode is governed by others and 
the time asymptotic response does depend on the 
coupling of Fourier modes- we must deal with the 
entire wave packet. 

Fortunately the system of Eqs. (|17tir5|) . with pe- 
riodic coefficients, is readily amenable to a sim- 
ple numerical solution for the densities N\ , N 2 
(and the electric field, E = N x - N 2 e ix ). Differen- 
tial equations with periodic coefficients (Mathieu 
equation being the simplest example) are known 
to display parametric instabilities in well defined 
regimes. Mathematica solutions demonstrate that 
model Eqs. (|17HT8|) are, indeed, paramctrically un- 
stable and acquire a hefty growth rate when b ~ 
uj/n (this behavior was predicted by the analytical 
dispersion relation). It turns out the time asymp- 
totic behavior, for cases of interest, does settle into 
a pattern that growth rate of the instability can 
be read off from the time evolution graph. We 
shall soon see that rotation (centrifugally) driven 
Langmuir waves can grow at a very fast rate. 

3. Discussion 

We shall now explore the detailed characteris- 
tic of the centrifugally driven Langmuir waves for 
the crab pulsar. Our eventual aim is to put to- 
gether a mechanism for delivering the rotational 
slow down energy to particle acceleration. It will 
happen in two major steps: (I) in the first stage, 
by means of the centrifugally induced parametric 
instability, the pulsar's rotation energy is trans- 
ferred to the electrostatic waves; (II) and in the 
second stage the amplified electrostatic modes un- 
dergoes Landau damping on high energy electrons 
causing efficient particle acceleration. 

3.1. Centrifugally induced Langmuir waves 

In order to make reasonable estimates for the 
instability growth rate, we recapture some details 
of our model. What we are studying is the insta- 
bility that arises in a highly relativistic electron- 
positron plasma in the pulsar magnctosphere. We 
have modeled the plasma in terms of two distinct 
streams distinguished by their 7 factors and the 
initial orbit phase (we have assumed the speed 
of both streams to be c). Let the two streams 
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have quite different 7s with 71 < 72 . The frequen- 
cies u)\ t 2 are calculated from the stream densities 
and 71, 2- Assuming that energy is uniformly dis- 
tributed throughout the magnetospheric plasma, 
we may approximate riiji w 71272 ~ ^fc7fc, where 
rib = n aj = VLB / (27rec£ 2 ) is the Goldreich-Julia n 
number density (jRuderman fc Sutherland! fl975h . 
The magnetic field (on the light cylinder zone) 
B = B st (R st /r ) 3 is related to B st 6.7 x 10 12 , 
the magnetic field close to the neutron star's sur- 
face. With these density relations, we conclude 
that one of the parameters for the numerical solu- 
tion, a = LO2 2 /lo\ 2 = 7i 4 /7| will be much smaller 
than unity. The most important determining pa- 
rameter of the system, however, is w = Wi/Cl, the 
ratio of the rotation frequency to the mode fre- 
quency (« wi); the parameter w will lie, typically, 
between 10 6 — 10 10 , the former near the light cylin- 
der and the latter near the pulsar surface. Speci- 
fying a, w, b, and the initial conditions is sufficient 
to seek a numerical solution. The parameter b is 
set by the wave number k and the phase difference 
0_ . From the general character of the differential 
equations, we expect strong parametric amplifica- 
tion when b > w. 

From our extensive numerical study, we pick 
up two representative Mathematica plots (Fig. [2] 
and Fig. [3]). In these, the solutions for ReNi, 
ImNi, ReN2, and ImN2 are plotted as functions 
of time. Time has been normalized to the inverse 
of the frequency uj\, the relativistic plasma fre- 
quency corresponding to the stream with lower 
7 and much higher density. The respective pa- 
rameters for the two figures are: For Fig. [21 



= b = 10 7 , a = .01, and for Fig. El 



10 s 



a = .01. Since the x axis denotes time measured 
in their respective plasma time (inverse of u>±)- the 
" absolute" time scales in the two figures arc differ- 
ent by a factor y/lO. For each case, we choose ini- 
tial conditions(t=0) for what we call the Sin solu- 
tion: iVi(O) = 0,AV(0) = l,iV 2 (0) = 0,JV 2 '(0) = 
a(cos b — i sin b); it is easy to verify that these con- 
stitute a consistent solution to the coupled differ- 
ential equations in the limit t — > 0. 

The principle results of the numerical investi- 
gations may be summarized as: 

a) For initial times t/w << l(not displayed in 
the figures), the streams are essentially uncoupled 
and each oscillates at its own plasma frequency; 
the frequencies are real. Since w is large, this 



stage can be reasonably long. For the two cases 
displayed, the period lasts « 3.10 4 plasma periods 
for case 1, and as 2.10 5 plasma periods for case 2. 
In each case the period is w 10~ 3 of the star rota- 
tion period. Compared to the star rotation time, 
the coupling and instability build-up time is very 
small for the displayed examples. 

b) As long as b <C w, even for large t/w, the cou- 
pling between the streams remains small. Conse- 
quently, the perturbations amplitudes remain lim- 
ited. 

c) As b approaches w, the two streams begin to 
interact after some initial lapse of time and all 4 
amplitudes begin to grow. Since the two streams 
are assumed to have different intrinsic plasma fre- 
quencies (a = 0J2 2 /0J1 2 = 7i 4 /7| << 1) the time 
dependence of (ReNi, ImNi), and (ReN2,ImN2) 
are quite different. It can be seen from the figures 
that although (ReNi, ImNi) oscillate (at a new 
slow frequency, the plasma oscillations are very 
fast and ever present-not resolved in the pictures) 
and grow, the amplitudes (ReN 2 , Im^) simply 
grow. Interestingly enough, the slow time evo- 
lution is exponential as predicted by the Floquet 
theory of differential equations with periodic coef- 
ficients. Thus the amplitude of entire perturbation 
(not just a given Fourier mode) exponentiates, and 
one can define a "global" growth rate. 

d) From the numerical solutions one can readily 
extract the the growth (exponentiation) time. For 
the solution in Fig. [2l it is « 2000 plasma times 
ss 2.10~ 4 P, while the growth rate for the solution 
in Fig. [3]is rj 6000 plasma times rj 6.10~ 5 P, where 
P is the pulsar period. Thus once the instability 
sets in, its growth rates are very high compared to 
the frequency of star rotation. 

e) For a given b and w, the growth rates go down 
as a 1 / 3 = (71 /72) 4 / 3 , precisely the scaling implied 
in the simple analytical estimate ([27]) . Given that 
the growth rates are quite high, one may conclude 
that even the streams with quite disparate 7 would 
produce a hefty instability. 

f) For a moderate and not too small a, the 
primary determinant of the instability gestation 
time, and growth rate is the relative value of 
b = 2cfc/Osin0_ and w = uji/Q,. Both toi and 
fl are determined by the pulsar phenomenology 
and are independent of the wave characteristics. 
Thus to bring b to the range of strong instability, 
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one has to choose fcsin0_ to be in the desirable 
range. In principle, then, for any finite sin</>_, 
one can find a range of k for strong growth. The 
centrifugal drive for Langmuir waves is, therefore, 
very robust and effective over a wide spectrum. 

g)Thc two stream model calculation, though far 
from perfect, is perfectly capable of capturing the 
essence of the pumping of Langmuir waves in the 
magnetospheric electron-positron plasma by the 
star rotation. A more detailed and comprehen- 
sive semi-kinetic model, that treats 7 and as 
phase space variables, is certainly recommended. 
However the current calculation creates a strong 
basis for the rapid build-up of rotationally pumped 
Langmuir waves (and their associated longitudinal 
electric fields) so that we can confidcntlly begin to 
explore the consequences of these electric fields for 
particle acceleration. 

3.2. Particle acceleration 

After having established that the rotational en- 
ergy of the pulsar can be rapidly converted into 
charge separation electrostatic fields, we will now 
address the next issue of driving fast particles 
with these electric fields. Looking for a mecha- 
nism to produce super- energetic particles, we will 
concentrate on further energizing the most ener- 
getic particles in the pulsar atmosphere- the pri- 
mary beam electrons that, to begin with, have 
rather high Lorentz facto rs as high as 7& ~ 10 7 
(jOsmanov fc Riegerl [20091) . 

The principal energy transfer (from waves to 
particles) channel explored in this work is the 
damping of Langmuir waves on electrons (Lan- 
dau damping). Landau damping is, primarily, a 
kinetic process and will need to be appropriately 
integrated into what has been, till now, an essen- 
tially fluid approach. Even at the cost of repeti- 
tion, it helps to emphasize that, in our model, 

1) The Langmuir waves are produced by the 
bulk e-p plasma characterized by the relatively 
lower 7 region of the distribution (Fig.l), 

2) while we want these waves to damp on the 
electron beam that appears as a small blip in the 
very high 7 end of the spectrum (Fig.l) 

The "producer" and the expected "consumer" 
particles, though sharing the same physical space, 
are very different in their numbers, and in per par- 
ticle intrinsic energy. 
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Y Y Y Y 

'p 't 'b ' 

Fig. 1. — The distribution function versus the 
Lorentz factor. As is clear from the plot, the func- 
tion consists of two parts: the first - wider region 
concerns the plasma component corresponding to 
the cascade processes of pair creation and the 
second one characterizes the primary Goldreich- 
Julian beam electrons. 



[Sin Solution, Nlr=S, Nli=Dot, N2r=BDash. N2i=SDash| 




Fig. 2. — The temporal behaviour of amplitudes 
oiReNx (S), ImNi (Dot), ReN 2 (BDash), ImN 2 
(SDash). 

(Sin Solution, Nlr=S. Nli=Dot, N2r=BDash, N2i=SDash| 













v\ 



Fig. 3. — The temporal behaviour of amplitudes 
of ReNx (S), ImNt (Dot), ReN 2 (BDash), ImN 2 
(SDash). 
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The study of electrostatic waves in ultrarel- 
ativistic plasm as started prior to the d iscov- 
ery of pulsars ( Siliu 1960t iTsitovich 1961 ). but 
gain ed considerable momentum after the discoy - 
ery (jLominadze et al. 1979; IVolokitin et al. I ll987ft . 



From this set of papers, one learns that (one di- 
mensional) electrostatic waves in relativistic e- 
p plasmas allow very simple dispersion in two 
limits: when the phase velocity is very high 
(v p h = uj r /k 3> c), and when the phase speed 
is very close to the speed of light, v p h ~ c. It 
is clearly the second case that is of interest to 
the energy transfer mechanism. In this limit, the 
real part of the dispersi on relation takes the form 
( Lominadze et al. I ll979l) 



oj = kc — P(k — fc ), 



(28) 



where /? = (7) / (27 s ) is a very small parameter () 
denotes momentum averaging), and fco is the limit- 
ing value of k for which the phase speed is exactly 
c. From Eq. ([28)) , we infer that the phase velocity 
will be less than the speed of light if k > fco; this 
condition defines the regime in which such a wave 
could strongly interact with physical particles. 

Going back to the unstable Langmuir waves, 
described in the previous section, we recall that 
the instability growth rates were found to be high- 
est when the parameter (argument of the Bessel 
function) (2cfc/f2) sin0_ is close to but moder- 
ately greater than [i res (the order of the Bessel 
function). The "resonance" condition, 

(2ck/Q) sin0_ « g/i r es = guj/Q (29) 

where g is number > 1, implies that the phase 
velocity of the unstable waves, u>/k = v p h ~ 
{2c/ g) sin0_, is in the correct regime as long as 
2 sin (/>_/<? is close to but less than unity. For g=l, 
for example, only the waves generated by streams 
whose phase difference, <f>- < 7r/6, will be able to 
accelerate particles; the highest growth pertain for 
waves for which </>_ is asymptotically close to ir/6 
and the phase velocity of such waves is asymptot- 
ically close to c. 

Thus the e-p bulk plasma, preferentially, gen- 
erates waves with phase velocity asymptotically 
approaching c. These waves will, preferentially, 
react with particles that move very close to the 
speed of light. In the pulsar magnctosphcre, the 
primary beam particles (the narrow blip on the 



high 7 end of the distribution function shown in 
Fig. [IJ, belong to the most relativistic class. The 
Langmuir waves, therefore, would, preferentially 
damp on them. Indeed, Eq. (f2"5]l tells us that 
closer the phase velocity to c, the smaller the value 
of /3 = (7) /(2j 3 ), and, consequently the larger the 
value of required (7) . 

To recapitulate the working of the Landau 
damping mechanism, let us place the phase veloc- 
ity of the wave in that part of the primary beam 
region in which the distribution function decreases 
with 7 (df /dp < or df/dj < 0). In this zone 
the number of electrons that satisfy u& < v p h is 
more than the number of electrons with Vb > v p h, 
Vb denoting the electron speed. Less energetic 
particles will gain energy from the waves, whereas 
more energetic electrons will loose energy, but 
since the number of particles with Vb < v p h is 
more than that of particles with Vb > v p h, the 
waves will damp with a net energy transfer to the 
beam electrons. 

As an aside, the parametric instability, dis- 
cussed here, may be different from the one con- 
jured in a typical laboratory setting. In the pulsar 
magnetospheric situation, the plasma frequency 
exceeds the rotation frequency by many orders 
of magnitude. The resonance is possible only for 
higher order harmonics and with phase velocity 
slightly less than the speed of light. 



From the relativistic plasma literature (jVolokitin et al 



19871 ). we can read off the estimated characteristic 
damping rate (jVolokitin et al.lll987l ) 



rvyu> p 
ni7i 



(30) 



where n, 7 and ui p = y / 47re 2 n/m are respectively 
the number density, the Lorentz factor and the 
plasma frequency of the specie on which the damp- 
ing occurs. 

If the instability growth rates, and the Landau 
damping rates were very disparate, there will be 
little effective energy transfer from waves to the 
particles. On the other hand, if the damping rates 
were far in excess of the growth rates, the waves 
will not grow much, again resulting in very little 
transfer from the star rotation to the waves. The 
most optimum scenario for an overall efficient en- 
ergy pumping/transfer system is realized when the 
instability growth and Landau damping rates are 
large and comparable, T ~ T^p. 
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Since it is the electron-positron plasma that 
supports and feeds the two-stream instability, the 
corresponding characteristic energy in the unsta- 
ble waves is of the order of the energy accumulated 
(transferred from star rotation) in the plasma 
components. Remembering that the number den- 
sity associated with the e-p-plasma far exceeds 
that of the beam component (see Fig.l), the en- 
ergy gained per beam electron could be extremely 
large if a substantial fraction of total "available" 
energy in plasma components were transferred via 
the electrostatic fields. 

The condition T ~ ^ld is satisfied, for in- 
stance, by the choice 72 = 4700, 71 pa 800 for 
the two streams. The resulting instability growth 
rate is of the order of 10 4 s _1 ; the corresponding 
growth period of 10~ 4 s is several orders of magni- 
tude less than the pulsar rotation period. The rate 
of energy transfer to the waves is extremely fast- 
the processes under consideration are extremely 
efficient. 

Let us probe further into the details of energy 
transfer from the star rotation into the plasma 
particles. In the local frame of reference the par- 
ticles are forced to slide along the magnetic field 
lines by means of the centrifugal force. Viewed in 
the laborator y frame, the reacti on force F reac pa 
2rocO^(r)- 3 ptogava et al.ll2003j) . driving the par- 
ticles, becomes infinite on the light cylinder sur- 
face. This is a natural result- to preserve rigid 
rotation, the particle velocity, in this region, must 
exactly equal the speed of light. Note also that 
the radial velocity tends to zero, because on the 
light cylinder the particles can only rotate with 
the linear velocity c. It is clear that the maximum 
energy gained from the rotator can be estimated 
as the work done by the reaction force. For all 
particles involved in the process, the energy gain 
may be estimated as (n^ 3> n 2 ) 



tides, ribSVe, 



W pa ni SVF reac Sr, 



(31) 



where SV is the volume in which the pumping 
takes place and 5r ~ c/T is the corresponding 
lengthscale. 

Since the aforementioned work done is trans- 
ferred to the beam electrons in the same volume, 
the energy of a given beam particle e can be ap- 
proximately calculated by equating W with the 
total energy gained by the accelerated beam par- 



niF reac 5r 
n b 



(32) 



One can readily show that, in the vicinity of the 
light cylinder surface (£ ~ 10~ 2 ), and for typical 
magnetospheric parameters, n\ pa 2.5 x 10 u cm~ 3 , 
Tib = n GJ ~2x 10 7 cm~ 3 , e, the energy acquired 
by a beam electron could get to be as very high- 
100s of TeVs or higher. 

Because the highly relativistic particles will, in- 
evitably, loose energy due to radiation, one must 
investigate how radiation will affect overall energy 
transfer. Could radiation, for example, put a strin- 
gent limit on the maximum energy acquired? 

For highly relativistic electrons and photons 
with ee p / l /(m 2 c 4 ) > 1 (e p h is photon energy), 
the inverse Compton mecha nism operates in the 
Klein-Nishina (KN) regime! Rvbicki &; Lightrnanl 
19791) . Energy emitted per particle per second 
is, then, given by the approximate expression 
Pkn c< TTr^m 2 c 5 n v h(eph) I I n (4ee/m 2 c 4 — 11/4) | 
( Blumenthal fc Goul (III 11)7(11) where r e is the classi- 
cal electron radius. The corresponding time scale 
of the process is tKN ~ c/Pkn- Using the ap- 
proximate photon number density, n p h(e P h) ~ 
L/4TrRf c ce p k (Ri c = c/Q is the the light cylin- 
der radius), and taking into account that for the 
Crab pulsar, L ~ 10 36 erg/s, e ph ~ lOOGcV, it is 
straightforward to show that if e = lOOTeV (for in- 
stance), the aforementioned timescale exceeds the 
typical instability time scale (ti n ~ 1/r) by many 
orders of magnitude. Thus Compton cooling is 
very slow compared to the wave energy transfer 
time. The KN time scale goes up with the par- 
ticle energy implying that at higher energies, the 
inverse Compton mechanism in the KN channel is 
tot slow to impose any constraints on the maxi- 
mum attainable electron energy. 

The next possible loss mechanism for relativis- 
tic particles, moving in magnetic field, is the syn- 
chrotron emission wi th the corresponding powe r 
P syn Pa 2e 2 w|7 2 /3c (|Rvbicki fc Lightmanl [l979h . 
where lob = eB/(mc) is the cyclotron frequency. 
One can readily show that the electrons leaving 
the gap with a 7 ~ 10 6 will undergo so effi- 
cient synchrotron cooling that the corresponding 
time scale, t syn ~ jmc 2 / P syn will be of the or- 
der of 10 _21 s. Thus, shortly after the particles 
leave the gap, they radiate their transverse mo- 
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mentum, and very soon transit to the ground Lan- 
dau state. After that, zipping only along the field 
lines, the electrons will reach the light cylinder 
zone in due course of time. It is precisely the re- 
gion, where the Langmuir waves, always propa- 
gating along the field lines, are excited. Therefore 
the wave interaction with particles will not cause 
pitch angle scattering, efficiently suppressing the 
synchrotron mechanism. Quasi linear diffusion, 
anot her possible source for imparting a pitch an- 



gle, (jLominadze et al. also does not oper- 



ate because the required condition for diffusion 
wb > 0)1,2, is violated for extremely energetic plas- 
mas (for plasmas with energy density exceeding 
that of magnetic field). The synchrotron mech- 
anism, therefore, is not expected to impose any 
constraints on the maximum attainable energies. 

How about the curvature radiation emitted by 
particles moving along the curved magnetic field 
lines? High energy electrons with energies in (100s 
of TeVs), i.e, with particle energy density exceed- 
ing magnetic energy density by many orders of 
magnitude, move, practically, in a straight line. 
This in turn excludes the curvature radiation from 
interfering with our energy transfer mechanism. 
As in previous cases, curvature radiation "does" 
not impose significant limits on e. 

The absolute upper bound on the energy that 
an electron may acquire, may be estimated as fol- 
lows: The total power budget, the slowdown lu- 
minosity, is 5.5 x 10 38 erg/s, about 40% of which is 
radiated away by the Crab nebula. If the remain- 
ing 60%), AL ~ 3.5 x 10 38 erg/s, were fully used 
up (via the mechanism proposed in this paper) to 
accelerate electrons to ultra high energies, a beam 
electron could emerge with an energy 



AL 



4nriRf c nGjc 



1- 



PeV 



V 



(33) 



where r\ < 1 is the fraction of particles involved 
in the process. Since the particles involved in the 
acceleration-radiation mechanisms, are from the 
region of open magnetic field lines, it is natural to 
assume that 77 <§C 1 ; this bunch of particles could 
be boosted up to multi PeV range. This estimate 
should be taken for what it is- an absolute upper 
bound. The most important message is that even 
if a very small percent of the total available energy 
could be transferred to fast electrons, the pulsar 
could be a source of ultra rclativistic electrons with 



energies in the range of 100s of Tevs. 

It is clear from the expression of F reac that the 
force increases asymptotically as we approach the 
light cylinder zone £ — > 0. However, the projected 
singularity at the light cylinder surface may be cir- 
cumvented by a recently conjectured mechanism: 
By means of a curvature drift instability, the mag- 
netic field lines twist' and, lagging behind the rota- 
tion, lead to the force-free dynamics of particles. 
The rigid rotation and the subsequent violation 



of the causality principle (jOsmanov et al. 2009) 



is, thus, prevented. As a result, the process will, 
then, stop due to the aforementioned instability. 

4. Summary 

1. We have constructed a simple but nontrivial 
theoretical framework to explore the possi- 
bility of particle acceleration driven by the 
rotational slow down of a pulsar. The first 
step in the proposed mechanism is the con- 
version of the rotation imparted centrifu- 
gal particle energy, via a "two stream in- 
stability", to electrostatic Langmuir waves 
in the electron - positron plasma residing in 
the pulsar magnctosphcre. Landau damping 
of these centrifugally excited electrostatic 
waves on the high energy primary beam elec- 
trons transiting through the same region, 
constitutes the second step. 

2. By carrying out a linear instability analysis, 
relevant to a non-autonomous system, it was 
shown that, for physically reasonable pa- 
rameters (pertinent for example to the Crab 
pulsar), the characteristic timescale for the 
growth of Langmuir waves is much less than 
the typical rotation timescales of particles. 
This means that the invoked instability is 
extremely efficient and can rapidly convert 
the star slow down energy into the electric 
field energy 

3. Landau damping of the unstable Langmuir 
waves on primary beam electrons that con- 
verts the electric field energy into particle 
energy is also shown to be equally rapid. 
The combination creates a very efficient 
"machine" that generates ultra high energy 
(multi TcV to PcVs) cosmic ray electrons; 
the star slow down luminosity being the ul- 
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timatc accelerator. 

The mechanism, developed in this paper, is 
rather general. The magnetospheres of active 
galactic nuclei (AGNs), for example, be just as 
suited to support centrifugally driven electrostatic 
waves. And these waves, could, equally efficiently, 
accelerate any charged leptons or hadrons. 
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